# generate draws from a standard normal distribution
y <- rnorm(n = 1000, mean = 0, sd = 1) |> round(3)
# show some draws
head(y)[1] 0.932 0.931 1.317 -0.380 -0.861 -0.827
[1] -0.0148
[1] 1.005713
2.5% 97.5%
-1.998225 1.977025
[1] -1.959964 1.959964